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Abstract 

In this paper, we extend our previous discussion of the Bose metal to the held 
tuned case. We point out that the recent observation of the metallic state as 
an intermediate phase between the superconductor and the insulator in the 
field tuned experiments on MoGe films is in perfect consistency with the Bose 
metal scenario. We establish a connection between general dissipation models 
and gauge field fluctuations and apply this to a discussion of scaling across 
the quantum phase boundaries of the Bose metallic state. Interestingly, we 
find that the Bose metal scenario implies a possible two parameter scaling 
for resistivity across the Bose metal-insulator transition, which is remarkably 
consistent with the MoGe data. Scaling at the superconductor-metal transi- 
tion is also proposed, and a phenomenolgical model for the metallic state is 
discussed. The effective action of the Bose metal state is described and its 
low energy excitation spectrum is found to be u> oc k 3 . 
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I. INTRODUCTION 



In a recent paper , we argued that a system of interacting Cooper pairs may form a 
gapless nonsuperfluid liquid, i.e. a metallic state, in two dimensions at T = 0. We called 
this state a Bose metal(BM). Although such an idea might seem rather counter- intuitive 
offhand since a system of delocalised bosons typically forms a superfluid(SF) state at T = 0, 
closer thinking suggests that it is not. A Bose metal is possible in the same sense that a spin 
liquid is possible in quantum antiferromagnets. In fact, we showed that it is just another 
variety of a spin liquid: an E 2 spin liquid. 

In amorphous superconducting films, particularly MoGe, metal-like states analogous to 
the ones in granular superconductors have been observed in the presence of a magnetic 
field. The fact that such a state exists as an intermediate phase between a superconduc- 
tor(SC) and an insulator (INS), as predicted in Ref. |J, is already a good hint that the Bose 
metal scenario might be applicable to these systems as well. In this paper, we generalize the 
arguments which were previously developed for the case of zero applied magnetic field, to 
the field tuned case. We show that the results of field tuned experiments in superconducting 
films, where such a metallic state has also been seen [fJ-[§], also fit in with our concept of a 
Bose metal. 

We first discuss (section II) the fact that gauge field fluctuations present in the vortex 
system in the JJA model play a key role in the formation of the BM state and show how 
generalized dissipation models can also similarly lead to an analogous BM state. We then 
discuss how the concept of the BM state may be generalized to the field tuned case. The 
observation of three phases in ref. f| is consistent with the Bose metal scenario. This leads 
to the suggestion of new scaling formulae for resistivity across the SC-BM and BM-INS 
transitions. These are presented in section III and are shown to work very well for the 
MoGe data. In section IV, we develop a phenomenological model for the metallic state 
based on the physics described above, which we find is in good agreement with the MoGe 
data. Finally in section V we discuss the effective action and the low energy excitation 
spectrum of the Bose metal state. Here, we also address the issue of why the BM state has 
not been observed in the previous analyses of the JJA model. 

II. THE JJA MODEL AND EFFECTS OF DISSIPATION. 

Dynamical gauge field fluctuations: In ref. 0, we considered a Josephson junction array 
(JJA) model with onsite and nearest neighbour repulsion for large Cooper pair fillings and 
zero external gate voltage (i.e. J2i (<%) = 0, where 5hi is the charge fluctuation operator 
at the site i |]I|). We showed using the duality transformation relations that this model 
maps onto a too-component quantum plasma of vortices^) and antivortices(^), i.e. a set 
of non-relativistic bosons, moving in a dynamically fluctuating gauge field A 11 . 

This picture of vortices moving in a fluctuating gauge field is a simple quantum mechan- 
ical extension of the results of classical phase fluctuations in a 2-d superconductor. In the 
classical regime, a 2-d superconductor (with phase fluctuations) maps onto a two compo- 
nent classical plasma (VV) undergoing screening by a static electric field (E = — VA ), as 
described by Kosterlitz and Thouless(KT) 0. In the quantum regime two new effects are 
important: 
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(a) the Bose statistics of the vortices and antivortices allows for the possibility of their su- 
perfluidity; and 

(b) because of quantum fluctuations, the electric field is no more static but dynamical, i.e. 
E = — Vv4 — (l/c s )dA/dr and is associated with a corresponding magnetic field B = V x A. 
Hence, expressed in terms of vortices, in the quantum regime, two processes are important: 
depairing of VV pairs (or, blowing up of vortex loops in (2+l)-d) and Bose condensation 
of depaired vortices and antivortices, corresponding to the loss of phase order and growth 
of charge order for the Cooper pairs, as discussed in ref. [|lj. The presence of dynamically 
fluctuating gauge field A causes the vortices and antivortices to pick up random Aharonov- 
Bohm(AB) phases exp(i J A.dl) which makes these two processes distinct from one another. 
If there were no fluctuating gauge field, then because of the quantum zero point motion, the 
vortices and antivortices would Bose condense into a phase coherent state as soon as they 
unbind, i.e. the film becomes insulating. Now because of the gauge field A, the situation 
here is different. If the gauge field fluctuations are small, the dephasing induced by random 
AB phases is weak and superfluidity of V and V is retained. On the other hand, if the gauge 
field fluctuations are very strong, then, because of the random AB phases, the unbound V 
and V fail to produce a phase coherent state, and the system is a non-superfluid liquid. 

The effects of these gauge field fluctuations may be seen very clearly in terms of a world 
line picture for the vortices(appendix A). We represent this as a bosonic system interacting 
with a fluctuating gauge field A. The partition function is Pl[7]||: 




where Sq (A) denotes the gauge part of the action and P is the permutation of the particles 
|[25|| . This type of partition function may be obtained, for example, by substituting j°(r, r) = 
J2iQiH r ~ r «( r )) an d j a ( r , T ) = J2i Qif? {r)5(r — rj(r)) in eqn.(15) in ref. 0], where qt is the 
vortex charge. The AB phase factor exp(z J A.dfi) appears explicitly here. When the gauge 
field fluctuations are weak, this factor is close to 1. So, this case is very similar to the 
Bose system without any gauge field, and the ground state is an entangled liquid, i.e. a 
superfluid (appendix A). Now, when the gauge field fluctuations become very large, the AB 
phase trapped by the bosons is of order 7r (modulo 27r), and the phase factor appears in 
Z B a with fluctuating signs. This would cause cancellation of several terms in Z B a when 
entangled configurations are present, implying that such configurations enter the partition 
function with a low weight and are correspondingly high energy states, just like the fermionic 
case (appendix A). So, in this case, the ground state is a disentangled liquid, i.e. a non- 
superfluid. Phase separation is not possible because of the long-range interactions present in 
the (original) vortex Bose system. Hence, as gauge field fluctuations increases, there should 
be a phase transition from the superfluid to a non-superfluid state. 

This discussion makes explicit why random AB phases causes disentanglement. The 
non-superfluid state of bosons, thus obtained due to dephasing by random AB phases, is 
a liquid rather than a solid where the particles are localised by fluctuations, because the 
gauge field fluctuates with time. Let us say that the vortices get localised by fluctuations. 
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Then, since they see a magnetic field that fluctuates with time, by Faraday's law, they 
see an induced electric field as well. This electric field will knock the localised particles 
out of their positions. Thus, the presence of a dynamically fluctuating gauge field in the 
quantum vortex system leads to the possibility of a non-SF liquid phase, and by duality |L0 



a consequent possibility of a metallic state for the Cooper pairs, the original players in the 
system. Because the vortices feel retardation effects due to gauge field fluctuations, they 
dissipate even in a pure system, thus leading to finite resitivity at T = 0. 

In summary, the effect of AB phases induced by the gauge field fluctuations is to disorder 
the vortices and to create a new variety of quantum liquid, some signatures of which we 
discuss in this paper. 

Effects of dissipation: The following four features about the Bose metal scenario emerge 
from the description of the JJA model above and in Ref. |L|: 

(i) when the Bose metal is observed, it exists as an intermediate phase between the SC and 
the INS; 

(ii) the observation of the BM phase is associated with two phase transitions: one from SC 
to BM and the other from BM to INS; 

(iii) the BM phase is dominated by dynamical gauge field(gf) fluctuations, as felt by the 
vortices; and 

(iv) there is a competition between gf fluctuations and quantum zero point motion of the 
vortices at the BM-INS phase boundary: in the INS phase, quantum zero point motion wins 
and the vortices are in the SF phase. 

It has been suggested by Mason and Kapitulnik that dissipation effects will quite 
generally help the formation of a metallic like state for vortices in the quantum regime. We 
show here that the effects of a generalized dissipation model on the quantum motion of the 
vortices also maps onto a set of non-relativistic bosons moving in a dynamically fluctuating 
gauge field, just like the JJA model. 

Consider a generic dissipation model (no static disorder): 

S= drd 2 r — | j(r) \ 2 + / drd 2 rd 2 r' p(r,r)V(r — r')p{r' , r) 



o In Jo 

+ \f drdr'd 2 rr(r, r)n a(3 (r -r',r- r')f(r, r') (1) 
2 Jo 

where j a {r,r) = ^iqirf{r)8{r — rj(r)), p(r, r) = Y^ili^ij — r i( T )) an d n = average vortex 
density. g« = vortex charge, and the partition function is Z = J Dri(r)e~ s . The first term 
is the quantum zero point motion term (also, known as the mass term); the second term is 
the long-range (logarithmic, in a pure system) interaction among the vortices and the last 
term is the dissipation term. Summation is implied in the dissipation term and a, (5 refer 
to the space coordinates. This action can be recast, by virtue of a Hubbard Stratanovich 
transformation on the associated partition function, as: 

S = J drd 2 r— \j(r) \ 2 + J dTd 2 rd 2 r'p(r,r)V(r — r')p(r',r) 

+ i f d 2 rdrj a {r, r)a a (r, r) + f drdr' 'd 2 rd 2 r' 'a a {r, T)K a/3 (r -r',r- r')a /3 (r / , r') (2) 



with K a p(oj ni k) = ri a p(u n ,k) (u n = Matsubara frequency). The dissipation kernel K a p 
here is model dependent and the field a M is introduced through the transformation. The 
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Maxwellian type of coupling j.a implies that represents a gauge field. For example, for 
the Caldeira Leggett heat bath K a/3 (uj n ,k) = 5 a/3 | u n \ /r), which is time-dependent. The 
dynamical nature of the gauge field follows from the fact that dissipation is associated with 
velocity dependent forces. 

Clearly, this mapping implies two things. Firstly, as the vortices move in the presence of 
dissipation, they trap a random AB phase exp(i J a.dl), and consequently the arguments of 
dephasing by random AB phases presented earlier in this section imply a disordered liquid 
phase for a dissipative system as well. For example, if we consider the case of Caldeirra- 
Leggett heat bath above, we see that: when dissipation 77 is much weaker than the quantum 
zero point motion, gauge field flucutations are weak, the randomness induced by the AB 
phases is very small and the delocalised vortices are in an entangled state, i.e. a superfluid. 
However, when dissipation is very strong, gauge field fluctuations are also correspondingly 
very strong, and the random AB phases can dephase the vortices into a disentangled, i.e. 
non-superfluid, state. The liquidity of this state follows from the dynamical nature of the 
gauge field. Because of this mapping, we shall quite often use gauge field fluctuations 
and dissipation interchangeably in this paper. Secondly, the above mapping means that 
a dissipative quantum vortex system (where vortices have been induced by an external 
magnetic field) has a very similar phase diagaram to that of the JJA model considered in 
Ref. except that now dislocations are present (see below). There will be an SC phase 
which consists of dislocation - antidislocation pairs, when quantum fluctuations are weak. As 
quantum fluctuations are increased (by tuning the field, for example), the pairs will unbind. 
The state this film would enter depends on the strength of dissipation. If dissipation is 
weak, the vortices will form a superfluid state and the transition is from SC to INS. On the 
other hand, when the dissipation is strong, the vortices would first enter an uncondensed 
liquid like state due to dephasing effect by the AB phases. By duality, this state is a non-SF 
liquid for the Cooper pairs (see first part of this section), and hence, it is fair to call this 
a Bose metal. The resistance is induced by the free dislocations in this case. As the field 
is increased further, quantum zero point motion overcomes dissipation (or, gf fluctuations), 
the vortices would Bose condense and the film will be insulating. Thus, the transition is 
SC-Metal-INS in this case. We explore some of the consequences of this scenario below. 

There are two key points which need to be kept track of while generalising the results for 
vortices from the zero field case to the field driven case in a dirty film, (a) In the presence of 
an external field, vortices enter an SC film in the form of an Abrikiosov lattice. True long 
range order is not possible in a 2d system at finite temperatures: dislocation-antidislocation 
pairs are created in the Abrikosov lattice. As temperature is increased, these pairs unbind 
(the analog of VV unbinding in the zero field case), and (Cooper pair) superconductivity 
is destroyed. (b)In the presence of static disorder, the Abrikosov lattice is converted into a 
glass, called a vortex glass(VG). Energy barriers between the various metastable states in 
this glassy phase are finite in 2d [p]-13|. One way to understand this is that dislocations are 



point like objects in 2-d, and that disorder screens long-range log interaction between the 
dislocations. Thus, the energy barriers to create dislocation-antidislocation pairs is finite in 
2-d, and the energy barrier to their motion is very small, particularly in the collective pinning 



regime ||13|| , in which we mostly focus on here. This means that as soon as the dislocation- 
antidislocation pairs are created, they will hop around and induce finite resistance at any 
finite temperature ~ exp(— €d/kT), where = energy barrier to creation of dislocation 
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pairs. This would imply that true superconductivity should set in at T = 0. However, at 
T = 0, tunneling processes resulting from quantum fluctuations due to quantum zero point 
motion, etc. can in principle be strong and destroy long-range order in the vortex glass 
phase, which is discussed in this paper. 



III. SCALING AT THE ZERO TEMPERATURE PHASE TRANSITIONS 

The presence of a Bose metallic phase in a SC film is associated with two phase transi- 
tions: one from SC to BM and another from BM to INS. Correspondingly there will be two 
scaling behaviours for the resistivity even for the field-tuned case, which we discuss below. 
1. SC-BM transition: The first transition is associated with the unbinding of (quantum) 
dislocation - antidislocation pairs (or, in a dirty system, when the free dislocation - an- 



tidislocation pairs come into existence |T3[)- The film enters a metallic state due to strong 
gauge field fluctuations. Finite resistance in the film is induced by free dislocations, which 
is proportional to the free dislocation density n$ |13|: R n ~ Rqu^iiv Rq = h/4e 2 is 



the quantum of resistance and fi v = vortex mobility, Udf scales as n<z/ ~ !/£+• £+ = SF 
correlation length that diverges across the SC-BM phase boundary with an exponent uq: 
£ + ~ (H — H c0 )~~ u °; H c0 is the critical field for SC-BM transition. Hence, on the metallic 
side, 

R D ~(H-H c0 ) 2 »°. (3) 

This scaling formula is quite different from traditional quantum SC to non-SC scaling R n = 
f{5/T 1 / uz ). Observation of this scaling in ref. H provides good evidence that the metallic 
phase is a Bose metal. A comment about the value of H cQ : since energy barriers to metastable 
states are finite in 2-d in the presence of disorder, one might think that a small amount of 
quantum fluctuations (zero point motion and dynamical gauge field fluctuations) would 
destabilise the VG order completely, i.e. H c0 = 0. But, at very low fields, the vortex 
system moves into the individual pinning regime from the collective pinning regime, and 



the energy barriers to the vortex motion due to pinning become very large []13" |. Hence, it 
is not inconceivable that, in this very low field regime, at T = 0, quantum fluctuations and 
disorder can conspire to produce a pinned vortex state without any free dislocations, and 
consequently true (Cooper pair) superconductivity. In other words, there is a possibility 
that H c0 > §. 

2. BM-INS transition: As the field is increased further, quantum zero point motion of the 
vortices (~ h 2 n v /m v ~ H) increases, and beyond a critical value H c , the zero point motion 
overtakes the gauge field fluctuations, the vortices form a SF phase and the film is insulating. 
Presence of free vortices during this transition motivates a scaling for the resistivity across 
the BM-INS phase boundary. So far, this phase boundary has been thought to be a SC-INS 
transition and the predicted scaling F3] is R n = f(5/T 1 ^ uz ). Although this scaling formula 



works at high temperatures, it fails at low temperatures 0. The two parameter scaling 
formula which we obtain below for the BM-INS transition, however, scales the data across 
the entire temperature range, both high and low, when the external field is in the critical 
regime. 

In ref. QXJ] , we mentioned that there is a jump in the vortex superfluid density at the 
BM-INS phase boundary. Naively, this would imply that this is either a first order or a 
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KT type transition. But, the transition is possibly more subtle than this because it is 
a phase transition from a (vortex) superfluid to a gapless nonsuperfluid phase. In a first 
order transition, there is no divergent length scale and the length scale is finite through the 
transition. But, the BM phase is gapless, i.e. characterised by a diverged length scale. The 
KT transition, which is a close precedent to this, on the other hand, is a transition from a 
superfluid to a gapped nonsuperfluid state. Also, the situation here is fundamentally different 
from KT because of the involvement of the gauge degrees of freedom. Thus, we suspect that 
this transition belongs to a very new universality class. In the absence of the knowledge of 
what exactly this universality class is, it is hard to write down the scaling formula across 
this BM-INS transition. Below, we suggest a way to scale the resistance across this metal 
insulator transition based on phenomenological considerations and leave the rest to future 
research. 

Since this is a second order quantum phase transition, there is a diverging correlation 
length £ with an exponent u at this transition, i.e. £ ~| 5 \~", and a frequency scale Q, which 
goes to zero with an exponent z, Q ~ where 5 = (H — H c ). As the energy dissipated 
scales as, (V 2 /R)t ~ Q (V = voltage drop, t represents the time), resistance R scales as 
R ~ V 2 /Q 2 . Since the dissipation is due to the vortices, the voltage induced by moving 
vortices is V — (h/2e)d8/dt, with dO/dt = 2nnfLv |TJ]], where v is the vortex velocity, L is 
the length over which the vortices move and n/ is the vortex density. To obtain the scaling 
of V, we note that the scaling of L is L ~ £ ~ Vt~ l l z . Since mv 2 ~ VL, the scaling of v is 
v ~ V"- We assume that the scaling oi rif is rif ~\ 8 \ a , where a is an exponent to be 
determined below, n/ is not the total vortex density H/<&q, but rather the critical fraction 
of the field induced vortices that participate in the dissipative process. Combining all these 
factors, we find that the scaling of R is R ~ 5 2a /VL 1+2 I Z . At any finite temperature, the 
divergence of £ is cutoff by temperature T, i.e. Q, ~ T. This also implies that the scaling 
function is always a function of 8/T 1 ' 1/z , i.e., 

RT 1+2/z /5 2a = f(5/T l/uz ) 

where / represents the scaling function. The right hand side comes from the fact that 
£ ~ 5~ u , ft ~ and Q ~ T, so that 5/T l l vz is the scaling variable. To obtain a, one 
needs to first note that the resistance saturates to finite values independent of temperature 
at low temperature; i.e. in the above eqn, we must have f(x) — > this limit. Plus, 

since this low temperature resistance of the film is non- critical through the metal insulator 
transition when H is tuned through H c as can be seen from eqn. (11) below, we should have 

2a = v(z + 2) (4) 

Thus, the scaling formula for the resistance is: 

•rpl/ 



R 



= f{5/T l ' vz ) (5) 



This is essentially a two parameter scaling formula, which is expected of a bosonic system 
|i~5fl . For MoGe films ||, z = 1 and v = 4/3, and hence, 2a = 4. The corresponding 



plot of eqn. (5) for these films is shown in fig. 1. The data collapse with this two parameter 
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scaling formula is quite remarkable. We think that this points out that there is a true metal- 
insulator transition at this critical field. Although resistance usually does not receive scaling 
under normal circumstances, we found at least one precedent to this in the literature, i.e. 
when a dangerously irrelevant variable is present WB . 



IV. PHENOMENOLOGY OF FIELD TUNED EXPERIMENTS 

In this part, we discuss the phenomenology of the metallic state in a dirty SC film, close 
to the metal-insulator transition, by setting up a quantum "pseudo-temperature" model || . 
This is based on the idea of competition between dissipation and quantum zero point motion 
of the vortices present near this transition point, presented so far. 



The energy barrier to create of dislocation pairs in a disordered vortex lattice is [13 
e d ~ ^ n {H c / H). At high temperatures, dislocations are created by thermal fluctuations, 
i.e. dislocation density is rid c ~ exp(— e^/ZcT), and the resistance induced by these is R n ~ 
RQ n d^v ~ e~ ed l kT . At low temperatures, quantum fluctuations take over, which we model in 
terms of a quantum (pseudo) temperature Tq: when Tq is small, the vortices being bosons 
are in a SF state, and when high, it is in a disordered/metallic state (real temperature 
remaining constant and small). Because whether the vortices are in a SF state or non-SF 
state is determined by the competition between the dissipation rj and quantum zero point 
motion, Tq should be a function of these two parameters. Since heating a system increases its 
entropy, and (a)enhanced quantum zero point motion (~ n v /m v ~ H) leads to superfluidity 
in a Bose system, and hence, reduced entropy, while (b)action of dynamical gf fluctuations 
(~ 77) increases (pseudo) entropy, we expect Tq, for not too small fields, to scale as 

Tq ~ (n v /m v r ~ HP (6) 

where a, (3 are constants > 0. (It can be any monotonic function of r]/H in general.) 
Currently, we do not know how to obtain a and f3. In what follows, we shall assume a = 
(3 = 1 and show that this quantum temperature model captures the general phenomenlogical 
features of the field-tuned expts in the aforementioned regime very well. 

Thus, expressed in terms of the quantum temperature model, as the field increases, the 
quantum temperature of the system decreases and when it falls below a temperature of 
the order of Kosterlitz Thouless temperature Tkt (since vortices here are 2-d bosons), the 
vortices would form a SF, i.e. the film would be insulating. Since T^t ~ h 2 n v /m v ~ H, 
the critical field H c above which the film is insulating is obtained when Tq ~ T^t, i-e. H c 
scales as (using eqn.(6)): 

ff c ~^~V^ (7) 

since r] scales inversely as R n . This makes sense physically, since as rj increases, gauge 
field fluctuations increase and it is harder to create a SF out of the vortices. This inverse 
dependence of H c on R n is broadly consistent with the trend observed in the MoGe films 
(see Table I in Ref. |7|). 

When Tq > Tkt, the vortices are in a non-SF state and the film displays a metallic 
response. From eqn.(6), the density of quantum dislocations, by analogy with the thermal 
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case, is then n d Q ~ exp(— ed/kTo). Thus, the resistance of this metal at a finite temperature 
is given by R a ~ R Q {n dC + n dQ )ii v , i.e. 

R n ~ ii! n [e- ed/feT + e ~ td/kT Q] (8) 

This simple additive estimate is a good first hand measure when both thermal activation 
and quantum tunneling processes are going on simultaneously JIB]] . Here we have used 
\x v ~ R n /RQ, which is true of the Bardeen-Stephen processes. This implies that when the 
real temperature T falls below Tq, the vortices enter the quantum regime. This determines 
the crossover temperature T cross at which the system moves from the classical to the quantum 
regime: 

T C ross ~T Q ~ n/H (9) 

Hence, as the field is increased, the crossover temperature decreases, as seen in the Mo- 
Ge experiments . This unusual crossover from the classical to quantum regime, which 
is very unlike that of a single particle system, results from the strong cooperative effects 
present in the vortex Bose system. 

Thus, from eqn.(8), in the quantum regime, the resistance saturates to 

Ra ~ R n e-^ kT * ~ R n exp(-C^^ln(^)) (10) 

where we have used dimensionless expressions for rj and H by rescaling them with Rq and 
H c2 respectively. C is a constant of order unity. Eqn.(10) has the same field dependence 



as that obtained in Ref. |H| if one expands the logarithm about the critical field H c , i.e. 

Ra ~ R n exp(-C^(^-^)), (11) 

It should be noted that the metallic resistance is non-critical across the metal-insulator 
transition. 

Thus, despite the simplicity of the quantum temperature model, it captures the basic 
phenomenology of the experiments very well. 

To summarise, we find several indications here, viz. (a)occurence of the metallic state as 
an intermediate phase between SC and INS phases, (b)scaling behaviour at the SC-Metal 
and Metal-INS transitions being in accordance with the Bose metal scenario predictions 
and (c)the fact that the phenomenology of the metallic state is as would be expected for 
a Bose metal state, to suspect that the metallic state observed in the magnetic field tuned 
experiments in the SC films is prrobably an adiabatic continuation of the Bose metal state 
presented in ref. [0. 



V. THE EFFECTIVE GAUGE FIELD ACTION OF THE BOSE METAL 

In this final section, we discuss the effective gauge action of the Bose metal based on 
our analysis of the JJA model in Ref. IJ. We obtain certain properties of this metal on the 
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basis of this action, including its low energy excitation spectrum. And then, we revisit the 
issue of why the Bose metallic phase has not been observed in the previous analyses of the 
JJA model or the Bose Hubbard (BH) model. 

Effective action: It is clear from the discussion in section II that a Bose metal is a liquid 
of uncondensed bosons (vortices) moving in a transverse gauge field. The effective action 
of this liquid can be obtained by integrating out the vortices in a one-loop approximation. 
This calculation has already been done in appendix E of ref. M: the transverse part of the 
gauge action can be read off directly from eqn.(E.3) there. Thus, the effective action of a 
BM is: 



1/3 



(12) 



where the coefficients a and b are finite and scale with free vortex density and may be read off 
from eqn. (E.3) in |l| and uo n = 2irnT(n =integer) are Matsubara frequencies. The spectrum 
of the longitudinal part of the gauge field follows from the fact that since the vortices 
and antivortices are free, they screen each other. Since the screening length scales as the 
density of free vortices rif and rif ~ l/£+, the above result is obtained. Physically speaking, 
the above spectrum of the transverse gauge field comes about, beacuse it is dynamically 



screened |22| and that there is no spontaneous symmetry breaking in this phase. Thus, the 
longitudinal modes are gapped but the transverse modes are gapless. 

Since the gauge field A 11 is seen by the vortices, one can use the duality transformation 
techniques []IU| , PT] on this gauge action to calculate the properties of the Cooper pairs (the 
original bosons) in the Bose metal phase. The superfluid density and the compressiblity of 
the charges, i.e. Cooper pairs, are [|l0| ,[2~lll : 

p c s = lim lim ^ \ /c2 =0 (13) 



q>_ J 

f^ow^o -fa^ + bq 2 c 



n c = lim lim — — — = — = (finite) (14) 



Here the superscript c denotes charges and we have used b = c 2 from eqn. (E.3), where c is 
the renormalised plasmon velocity. These results are consistent with the fact that a Bose 
metal is a non-superfluid compressible liquid. 

The specific heat of the BM can be calculated from eqn. (12). This has already been 
done in ref. [TJ. The gapless transverse modes make the most important contribution to the 
specific heat: the low temperature specific heat of the BM goes as C ~ T 2 / 3 . This T 2 / 3 
power law specific heat, in turn, directly implies that the low energy excitation spectrum of 
the original Cooper pair system in this liquid state goes as: 

u oc k 3 (15) 

This unconventional excitation spectrum results from the strong correlation effects present 
in the system and is consistent with the fact that Landau critical velocity should be zero in 
a nonsuperfluid Bose system. The density of states is correspondingly divergent: N(u) ~ 
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cj -1 / 3 . Obviously, the BM state cannotbe continued to a Fermi liquid or mapped onto a set 
of free particles. We shall argue below that this liquid is gapless beacuse it is an E 2 spin 
liquid, rather than an SU{2) spin liquid. 

In ref. 0, we pointed out that the phase ((pi) and charge fluctuation operators (Srii) 
entering the JJA model are the generators of E 2 , the Euclidean group in two dimensions; 
E 2 being a group contraction of SO(3). In this sense, the JJA model is an E 2 "spin" model 
and the BM phase being the disordered phase of this model should be regarded as an E 2 
spin liquid. An extremely important point is that, whereas an SU(2) spin liquid is usually 
gapped, an E 2 spin liquid is necessarily gapless. This follows from the distinction that 



whereas an SU(2) spin model maps onto bosons in a constant magnetic field [23], an E 2 spin 
model maps onto bosons in a fluctuating magnetic field. Because of the constant magnetic 
field, the SU(2) spin liquid maps onto a Quantum Hall liquid p3 |. This means that the 



effective gauge action has a Chern Simmons term and the excitations are correspondingly 
gapped. On the other hand, the effective gauge action of the E 2 spin liquid is of Maxwellian 
type, i.e. without any Chern Simmons term: 



S([X\) = J dud 2 qf(u,q)AZ^_ 



Now, if the excitation spectrum is gapped, i.e. lim w ^ f(u, q) = constant, then this means 
that there is spontaneous symmetry breaking and the z-component of the E 2 spins are 
ordered, i.e. the system is charge-ordered plUlfl- Hence, the E 2 spin liquid is necessarily 
gapless. For the JJA model we have been discussing so far, this feature can be seen directly 
by focussing on the transverse part of the gauge action in eqn.(12). 

It is hard to rule out the possibility of excitations with fractional quantum numbers 
H in this liquid. However, since there is no explicit time reversal symmetry breaking in 
the effective action, we think that such excitations would always exist in pairs being held 
together strongly by gauge-like forces. 

JJA model: what's missing? In this last part of the paper, we revisit the issue of why the 
Bose metallic phase has not been observed in the previous analyses of the JJA model or the 
Bose Hubbard (BH) model. 

The key issue which distinguishes our analysis from the others is that our calculation has 
been done in the limit when the average bosonic filling per site is very large. In this limit, 
the BH model maps onto a JJA model. A lot of current calculations on the JJA model are 
strongly influenced by the results available for the BH model at low fillings. However, as we 
stressed in our previous paper, the algebraic properties of the BH model in these two limits 
are very different: in the limit of low fillings, the model is very close to that of the SU(2) 
spin model. On the other hand, in the limit of large fillings, it is close to that of the E 2 spin 
model. 

It has been recently pointed out that SU (2) and the E 2 algebras are actually connected 
by a singular transformation [[2(J so that the results obtained in one limit cannot be continued 
adiabatically to those in the other limit, i.e. without crossing a phase boundary. This also 
lends extra support to our argument above that the properties of the spin liquids supported 
by the two spin models are distinctly different. Hence, it is very important to maintain this 
distinction between SU(2) and E 2 while performing calculations on the JJA model. The 
limit of large fillings is good for the SC films, whereas low filling limit is good for helium 



films 121 
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As a second point, a lot of the calculations on JJA model have been done for the case 
when an external gate voltage V g is present. We have worked on the case V g = 0. This 
situation is close to that of the real granular SC films. And, finally, our calculations have 
been done in the presence of nearest neighbour interaction V\. The presence of V\ brings 
about non-trivial effects in the JJA model, a quite well explored instance of which is the 
onset of supersolidity when rescaled gate voltage is close to half integer. More details on 
these differences may be found in Ref. [[2~T|]. 

There are two previous calculations on the JJA model for V g = and V\ 7^ 0: mean 
field calculations and numerical simulations |24|| . The reasons for non-observation of the 
BM phase in these are as follows: (a) Mean field: A BM phase benefits partly both from the 
kinetic energy and the potential energy terms. But, in a mean field calculation, as soon as 
superfluidity is destroyed, the kinetic energy term is completely suppressed, (b) Numerical 
Simulations: The problem with this calculation is that the authors investigated only the 
superfluid phase in the V — V\ plane at V g = 0. Only superfluid density and structure factor 
Sn,* were measured here, which do not distinguish between a Mott insulator phase and a 
BM phase: both of these quantities are zero in these states. More simulations addressing 
this issue might be helpful to crosscheck the existence of the BM phase and its properties. 
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APPENDIX A: WORLD LINE PICTURE FOR BOSONS AND FERMIONS 

Here we enlist the actions for bosons and fermions in the world line picture, which 
might help the reader in understanding the world line picture discussion of the gauge field 
fluctuations in section II. 

First, consider a system of N interacting bosons in the world line picture, without any 
gauge field. Partition function of this system is 



z B = -^Y,f (- f ^ 

N\ V Ari(/3)=r Pi (0)> V \ JO 



i i^j 



where P is the permutation of the particles ||25|| . At low temperatures, the world lines get 
entangled because of quantum zero point motion effects, which implies a finite superfluid 
density for the bosons 



Now consider an equivalent system of fermions [26]. Here, the partition function is 



JttP 



-Te 



dr 



i i^j 
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the extra phase factor coming from the anticommutativity of the fermions. Because of 
this phase factor, it is readily clear that entangled configurations enter with random signs 
and cause cancellation of several terms. Thus, since the entangled configurations enter the 
partition function with low weight, this means that these configurations are high energy 
configurations and disentanglement is favoured for a fermionic system. And, hence the 
ground state of a Fermi system is a non-superfluid. 

The case of bosons with gauge field fluctuations discussed in section II is intermediate 
between these two cases. 
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FIG. 1. Scaling collpase of the MoGe data from Ref.2 across the metal insulator transition 
on a log-linear plot using eqn.(5): 5 = H — H c ; H c = 1.22T is the critical field at which the 
metal insulator transition occurs. The data shown covers the entire low temperature range of 
measurement: ~ 0.02 — 0.2K. 
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